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I.  Introduction  and  Program  Objective 

Current  methods  of  lifetime  assessment  leave  much  to  be  desired.  Typically,  the  expected  life 
of  a  full-scale  component  exposed  to  a  complex  environment  is  based  upon  empirical 
interpretations  of  measurements  performed  on  microscopic  samples  in  controlled  laboratory 
conditions.  Extrapolation  to  the  service  component  is  accomplished  by  scaling  laws  which,  if 
used  at  all,  are  empirical;  little  or  no  attention  is  paid  to  synergistic  interactions  between  the 
different  components  of  the  real  environment.  With  the  increasingly  hostile  conditions  which 
must  be  faced  in  modem  aerospace  applications,  improvement  in  lifetime  estimation  is  mandated 
by  both  cost  and  safety  considerations. 

This  program  aims  at  improving  current  methods  of  lifetime  assessment  by  building  in  the 
characteristics  of  the  micro-mechanisms  known  to  be  responsible  for  damage  and  failure.  The 
broad  approach  entails  the  integration  and,  where  necessary,  augmentation  of  the  micro-scale 
research  results  currently  available  in  the  literature  into  a  macro-scale  model  with  predictive 
capability. 

In  more  detail,  the  program  will  develop  a  set  of  hierarchically  structured  models  at  different 
length  scales,  from  atomic  to  macroscopic,  at  each  level  taking  as  parametric  input  the  results 
of  the  model  at  the  next  smaller  scale.  In  this  way  the  known  microscopic  properties  can  be 
transported  by  systematic  procedures  to  the  unknown  macro-scale  region.  It  may  not  be  possible 
to  eliminate  empiricism  completely,  because  some  of  the  quantities  involved  cannot  yet  be 
estimated  to  the  required  degree  of  precision.  In  this  case  the  aim  will  be  at  least  to  eliminate 
functional  empiricism.  Restriction  of  empiricism  to  the  choice  of  parameters  to  be  input  to 
known  functional  forms  permits  some  confidence  in  extrapolation  procedures  and  has  the 
advantage  that  the  models  can  readily  be  updated  as  better  estimates  of  the  parameters  become 
available. 


II.  Program  Organization 

The  program  has  been  organized  into  specific  tasks  and  subtasks  as  follows. 

Task  100.  Lifetimes  of  metallic  dispersed-phase  composites 

Most  service  materials  fall  into  the  category  of  dispersion-hardened  metallic  composites.  This 
task  will  consider  the  problem  of  dispersion  hardened  materials  in  general,  but  with  two  specific 
materials,  NiAl  and  MoSi2/SiC  in  mind. 

Task  110.  Identification  and  modelling  of  micromechanisms 

The  purpose  of  th:s  task  is  to  determine  what  micromechanisms  are  operative  in  the  high- 
temperature  deformation  of  dispersion-hardened  materials.  In  the  general  case  this  will  be  done 
by  a  literature  search.  For  specific  materials,  the  micromechanisms  will  be  determined  from  the 
experimental  program  at  NRL.  Once  identified,  each  of  these  micromechanisms  will  be 
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modelled,  in  order  to  determine  what  are  the  critical  parameters  which  determine  its  effect  on 
plastic  flow  and  values  for  these  parameters.  Also  to  be  determined  is  whether  the  modelled 
critical  values  are  dependent  on  quantities  which  must  be  obtained  from  a  smaller  scale  model. 

Task  111.  Equiaxed  dispersoids 

This  task  will  consider  dispersions  of  the  type  encountered  in  NiAl-like  materials.  That  is,  the 
dispersoids  are  considered  to  be  small  compared  to  the  grain  size.  The  term  ’equiaxed’  is  used 
because  the  particles  are  roughly  of  the  same  size  in  all  three  dimensions.  However,  this  is  not 
a  requirement  for  this  task.  Rather,  it  is  necessary  that  the  particles  not  be  too  large  in  the 
dimension  normal  to  the  slip  plane,  so  that  they  can  be  surmounted  with  relative  ease  by  cross¬ 
slip  and/or  climb  without  the  generation  of  appreciable  back-stress. 

Task  112.  Anisotropic  dispersoids 

This  task  covers  the  case  of  dispersoids  which  are  elongated  in  the  direction  normal  to  the  slip 
plane.  An  example  is  SiC  fibers  in  MoSi2.  In  this  case,  plastic  flow  around  the  dispersoids  takes 
place  by  a  combination  of  glide  and  climb,  but  is  a  protracted  process  during  which  large 
stresses  acting  in  opposition  to  the  applied  load  are  developed. 

Task  113.  Grain  boundary  effects 

This  task  will  examine  the  role  of  grain  boundary  processes  in  high-temperature  deformation. 

Task  120.  Macroscopic  stochastic  mode!  for  creep 

In  real  materials  it  is  likely  that  more  than  one  mechanism  will  be  operative,  either  in  parallel 
or  in  series.  The  information  gained  in  task  110  is  not  sufficient  to  describe  this  situation.  Once 
the  critical  parameters  for  individual  mechanisms  have  been  determined,  it  is  necessary  to 
combine  them  in  a  macroscale  stochastic  model.  This  will  be  done  by  determining  critical 
stresses  and  activation  enthalpies  as  a  function  of  local  geometry  and  using  these  values  in  a 
finite-temperature  simulation  of  creep  through  a  random  array  of  dispersoids.  Careful  attention 
must  be  paid  to  possible  interactions  between  mechanisms. 

Task  130.  Extension  to  cyclic  deformation 

The  final  step  in  task  100  is  to  extend  the  results  to  the  case  of  cyclic  deformation.  Irreversibility 
is  an  intrinsic  feature  of  the  model  in  task  120.  However,  it  is  likely  that  other,  as  yet 
unrecognized,  characteristics  of  cycled  deformation  will  have  to  be  considered. 

Task  200.  Lifetimes  of  piezoelectric  ferroelectrics 

Failure  in  cyclic  loading  of  sensors  and  actuators  formed  from  lead  zirconate  titanate  (PZT)  is 
a  continuing  problem.  PZT  is  a  ceramic  and  therefore  differs  from  the  materials  considered  in 
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task  100  in  that  plastic  deformation  is  not  involved.  This  task  will  examine,  modelling  as 
necessary,  the  operation  of  PZT  devices,  in  order  to  determine  the  factors  governing  lifetime 
limitation. 

Task  300.  Reporting 

Running  concurrently  with  tasks  100  and  200,  this  task  will  inform  the  Navy  Program  Manager 
and  Contracting  Officer  of  the  technical  and  fiscal  status  of  the  program  through  R&D  status 
reports. 
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III.  Technical  Progress 
Task  111.  Equiaxed  dispersoids 
a.  Theory 

Hardening  of  metals  by  means  of  second-phase  or  extrinsic  particles,  commonly  known  as 
dispersion  hardening,  has  been  a  standard  route  to  improvement  of  the  mechanical  properties  of 
metals  for  many  years.  There  are  several  distinct  mechanisms  for  the  hardening,  all  of  which 
involve  retardation  of  plastic  flow. 

•Orowan  looping,  or  the  impedance  to  flow  offered  by  impenetrable  obstacles. 
•Stacking  fault  hardening,  in  which  the  stacking  fault  energy  inside  the  obstacle 
is  different  from  that  in  the  matrix. 

•Order  hardening,  for  which  the  crystal  structure  within  the  obstacle  forms  an 
ordered  superlattice,  so  that  the  low-energy  deformation  route  requires  paired 
dislocations. 

•Chemical  hardening,  which  results  from  the  surface  energy  generated  when  a 
precipitate  is  intersected. 

•Modulus  hardening,  due  to  a  mismatch  of  elastic  constants  between  dispersoid 
and  matrix. 

•Coherence  hardening,  caused  by  the  back-stress  from  inclusions  which  do  not 
fit  into  the  holes  in  the  matrix  without  elastic  deformation. 

For  each  of  the  mechanisms  listed  above,  there  exists  an  approximate  theory,  usually  supportable 
with  a  single  set  of  experimental  measurements.  The  experiments  themselves  are  difficult, 
because  the  mechanisms  are  not  separable  and  the  range  of  variables  which  can  be  accessed  is 
severely  limited.  This  is  not  at  all  a  satisfactory  foundation  on  which  to  build  a  global 
deformation  model. 

Our  unified  approach,  based  on  extensive  computer  modelling  of  individual  mechanisms,  has 
been  to  recognize  that  the  problem  is  essentially  geometric,  the  primary  factors  being  the  size 

D  and  spacing  L  of  the  dispersoids,  as 
illustrated  in  Figure  1.  The  detailed  nature  of 
the  particular  mechanisms  is  of  no  practical 
consequence,  and  can  be  accommodated  by 
assigning  an  intrinsic  strength,  or  impedance, 
to  each  obstacle.  This  intrinsic  impedance 
can  also  be  expressed  in  geometric  terms  by 
the  breaking  angle  6  (see  Figure  1).  The 
strength  of  any  dispersed-phase  composite 
then  depends  only  upon  the  statistical 
distributions  of  these  three  geometric 
parameters. 


Figure  1  Definition  of  the  geometric  parameters 
D,  L  and  $. 
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Work  to  date  has  concentrated  on  regular  dispersions,  that  is,  those  for  which  the  three 
parameters  L,  D  and  6  are  single-valued  throughout  the  material.  Even  this  limited  theoretical 
model  gives  a  better  description  of  the  available  experimental  work  than  any  earlier  treatment. 
Specifically,  we  have  reprocessed  or  performed  computations  of  the  strength  of  regular  arrays 
for  the  cases  of  coherence  hardening1,  modulus  hardening2,  stacking  fault  hardening3  and  order 
hardening3,  in  both  particle  intersection  and  Orowan  looping  limits.  We  find  that  in  all  cases  the 
strength  is  described  well  by  the  expression 


a  =  — sinG 
271 AL 


where  X  is  the  harmonic  mean  of  D  and  L,  that  is 


1=1 
X  ~  D 


(1) 


(2) 


and  A  is  a  constant  which  is  unity  for  screw  dislocations  and  (l-i')  for  edge  dislocations.  Use 
of  the  harmonic  mean  for  X,  which  is  the  characteristic  length  describing  the  elastic  energy  of 
the  dislocation  configuration,  reflects  the  dominant  influence  on  the  line  energy  of  the  dispersoid 
and  interdispersoid  dipoles  (formed  by  the  arms  XX  and  X’X,  respectively,  in  Figure  1). 

It  is  important  to  note  that  we  find  equation  (1)  to  be  valid  for  all  applied  stresses,  even  those 
less  than  the  critical  values.  This  means  that  (1)  applies  to  the  creep  (unidirectional  and  cyclic) 
conditions  important  to  lifetime  studies. 

b.  Comparison  with  Experiment 

There  is  no  question  but  that  the  fundamental  hardening  stems  from  an  Orowan-like  impedance 
and  therefore  depends,  to  first  order,  inversely  on  the  dispersoid  separation  L.  The  classical 
Orowan  model,  based  on  elementary  line  tension  theory,  predicts  the  dependence 


Figure  2  An  Orowan  plot  of  strength  for 
dispersion  hardened  materials. 


o  =  —  sin0  (3) 

L 

with  the  breaking  angle  factor  sin0  =  1. 
However,  a  simple  plot  of  the  hardening  A<r 
against  1/L  (Figure  2)  reveals  that  this 
approximation  is  grossly  inadequate.  The 
experimental  data  shown  in  Figure  2  consist 
of  79  measurements  from  eight  distinct 
sources4'11,  encompassing  all  of  the 
hardening  mechanisms  detailed  above.  The 
solid  line  indicates  the  theoretical 


6 


dependence  (3).  In  principle,  a  knowledge  of  the  true  breaking  angles  might  improve  the 
agreement,  but  none  of  the  experimental  work  includes  this  information,  which  would  require 
in-situ  TEM  measurements.  Suitable  values  for  0  could  bring  all  values  below  the  line  in  Figure 
2  into  agreement  with  (3),  but  not  the  many  points  which  lie  above  the  line. 

To  compare  the  data  with  our  unified  approach,  equations  (2)  and  (3),  we  use  the  correction  for 
coherent  precipitates  introduced  by  Duesbery  and  Sadananda1,  which  defines  the  effective 
diameter  D*  of  a  coherent  obstacle  by  the  relation 


(El)*  =  0.25Le  (4) 

D 

where  e  is  the  misfit  factor.  There  remains  the  problem  of  the  unknown  breaking  angle.  The 
approach  we  have  taken  for  this  is  to  use  the 
experimental  data  to  derive  the  distribution 
N(0)  required  for  the  experimental  points  to 
fit  the  theoretical  prediction.  The  results  for 
the  data  plot  are  shown  in  Figure  3.  Although 
the  fitting  of  the  breaking  angle  guarantees 
that  most  of  the  data  points  will  fall  exactly 
on  the  theoretical  line,  this  is  not  entirely  a 
trivial  exercise.  In  Figure  3  only  two  points 
lie  above  the  theoretical  line,  and  these  only 
just  so.  This  contrasts  favorably  with  the  case 
shown  in  Figure  2. 


The  distribution  N(0)  is  shown  in  Figure  4.  It 
can  be  seen  that  the  required  breaking  angles 
fall  mostly  in  the  range  of  10°  -  50°,  much  smaller  than  the  90°  usually  assumed,  but  consistent 

with  the  angles  observed  in  computer  modelling 
of  coherent2  and/or  penetrable2,3  dispersoids. 

We  cannot  state  unequivocally  that  the 
experimental  results  prove  that  our  unified 
approach  is  valid.  This  is  simply  because  the 
experimental  measurements  do  not  contain 
enough  information  for  this  purpose.  However, 
we  can  state  with  confidence  that  all  of  the 
experimental  work  is  consistent  with  the  unified 
approach. 


Figure  3  Plot  of  strength  following  the  unified 


approach  of  equation  (2). 


Figure  4  The  distribution  of  breaking  angles. 
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c.  Implications 

★  The  importance  of  the  unified  description  (2)  is  that  it  applies  to  all  passing  mechanisms  (with 
the  possible  exception  of  cross-slip  and  climb),  for  all  dispersed  phase  composites  in  which 
classical  dispersion  hardening  dominates.  The  utility  of  the  approach  is  that  it  permits  the  ready 
construction  of  a  global  geometric  model  in  which  the  dispersion  is  characterized  by  the 
characteristic  length  X  and  the  details  of  the  dislocation-particle  interactions  are  subsumed  into 
the  breaking  angle  distribution  N(0).  This  permits  the  ready  calculation  of  activation  enthalpies 
for  plastic  flow. 

★  The  unified  description  applies  not  only  to  ultimate  strength  estimation,  but  also  to 
deformation  by  unidirectional  and  cyclic  creep. 

d.  Plans 

★  Develop  this  unified  approach  into  a  stochastic,  high-temperature  creep  model  in  two  steps: 

1.  Creep  models  for  specific  mechanisms  (as  defined  by  the  breaking  angle)  with  both  regular 
and  random  arrays  of  dispersoids,  in  order  to  determine  stress  exponents  for  comparison  with 
and  analysis  of  experiments. 

2.  A  combined  model  for  creep  through  a  random  array  of  dispersoids  with  a  spectrum  of 
strengths. 

★  Hence  compute  lifetimes  in  unidirectional  and/or  cyclic  deformation  as  estimates  of  the  time 
required  to  reach  a  critical  level  of  damage. 
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Task  112.  Anisotropic  dispersoids 

An  analytic  model  of  the  hardening  due  to  elongated  fibers,  or  whiskers,  has  been  developed. 
a.  Theory 

We  assume  slender  cylindrical  whiskers  of  length  21  and  radius  r  to  lie  parallel  to  the  axis  of 
deformation  and  not  to  take  part  in  that  deformation.  The  first  assumption  does  not  impact  the 
conclusions  significantly.  We  further  suppose  that  the  material  in  which  the  whiskers  are 
imbedded,  strains  uniformly  at  a  rate  e.  This  deformation  is  achieved  by  processes  which  differ 
in  detail  with  distance  from  the  interface  of  the  non-deforming  cylinders.  At  remote  positions 
deformation  is  achieved  simply  by  dislocation  glide.  In  the  neighborhood  of  a  cylinder  and 
specifically  at  distances  from  the  surface  of  the  order  of  interatomic  spacings,  glide  is  inhibited 
and  is  replaced  by  ’glimb’  (a  combination  of  glide  and  climb,  possibly  alternating)  while 
individual  dislocations  form  elliptical  loops  which  circumnavigate  the  perimeter  of  a  particular 
cylinder.  Each  such  loop  retains  its  own  characteristic  Burgers  vector,  but  the  mean  direction 
of  all  such  Burgers  vectors  is  parallel  to  the  direction  of  straining.  Then  referring  our 
considerations  to  the  steady  state  in  which  incremental  strain  is  everywhere  continuous  we 
require  that  the  net  forces  and  motion  of  the  glimbing  dislocations  be  parallel  to  the  interface 
and  that  the  number  of  dislocation  of  magnitude  b  which  glimb  past  a  point  distant  x  from  the 
mid  point  of  a  whisker  be 


|  -  jwk*)  (5) 

where  v(x)  is  the  dislocation  velocity  and  f(x)  the  density  (that  is  the  number  of  dislocations  of 
magnitude  b  per  unit  length)  at  the  said  position  x.  The  force  acting  per  unit  length  (L)  of  any 
dislocation  is 


F  =  obL  (6) 

when  the  stress  o  and  the  Burgers  vector,  of  magnitude  b,  are  so  disposed  that  the  force  on  the 
dislocation  is  in  the  direction  of  the  Burgers  vector.  Here,  as  stated  above,  the  mean  direction 
of  the  Burgers  vectors  of  all  dislocation  loops  is  parallel  to  the  axis  of  strain  and  we  shall  in  the 
first  instance  suppose  that  this  direction  is  parallel  to  that  of  the  axes  of  the  cylinders.  Taking 
advantage  of  the  fact  that  this  model  requires  that  these  dislocations  move  parallel  to  this  same 
direction  we  replace,  fo  the  purposes  of  calculation  of  stresses,  the  actual  arrangement  of 
elliptical  dislocation  loops  with  a  continuum  of  circular  loops,  having  displacement  vectors 
parallel  to  the  mean  of  their  directions  and  this  to  the  axes  of  the  cylinders.  We  further  require 
that  the  glimb  of  these  dislocations  involve  the  emission  and  absorption  of  vacancies  and  thus 
the  processes  of  diffusion.  To  the  extent  that  the  velocity  resultant  on  the  application  of  a  force 
may  be  regarded  as  the  consequence  of  an  Einstein  drift  it  is  given  by 
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DF  _  Dba_ 
kT  "  kT 


(7) 


where  D  is  the  coefficient  of  self  diffusion  of  the  matrix,  k  is  Boltzmann’s  constant  and  T  the 
temperature.  Since  the  dislocation  motion  is  restricted  by  the  need  for  climb  it  can  be.  argued 
that  F  should  be  identified  as  the  force  available  for  that  motion.  However,  it  would  seem  likely 
that  the  actual  motion  would  be  clearly  separable  into  climb  and  glide  over  only  very  short 
lengths  (i.e.  a  few  Burgers  vector  distances)  with  the  result  that  the  force  available  would  be  that 
for  the  mean  motion.  The  stresses  developed  on  the  glimbing  dislocations  are  given  by 

°M  =  -g-J'MOM*  (*> 

2it(l -v)J  -i 

Because  the  dislocations  which  surround  the  cylinder  are  taken  to  be  circular  loops,  with  Burgers 
vectors  normal  to  the  planes  of  the  loops  we  see  that  the  stress  relevant  in  the  calculation  of  the 
force  parallel  to  the  axis  of  the  cylinders  is  that  for  shear.  The  stress  from  such  a  loop  has  been 
calculated  by  Kroupa12.  However,  one  may  as  first  approximation  avoid  the  usage  of  the 
complicated  representation  of  the  stresses  by  recognizing  that  for  small  values  of  t-x,  that  is 
when  the  factor  h  =  (t  -x)  /r  <  <  1,  that 

G(tj)  =  —  (9) 

t-x 


and  decreases  very  rapidly  as  h  increases  outside  this  range.  Then,  as  a  first  step  towards  a 
simple  and  reasonably  accurate  approximation  we  write  (8)  as 


o(x)  = 


\lb  far*x  fit)  ^ 
27t(l-v)J-<w*Jf  t~X 


(10) 


where  a  is  a  constant  of  order  unity.  Expanding  f(t)  in  a  Taylor  series  about  the  point  x  this 
becomes 


o»  — t b—rx  *  ...* 

27T(1  -vy-ar*x  t-X  dx  2  dx2 


(11) 


and  we  have 

OCX)  =  -  ^  V  (12) 

75(1  -V)  dx 


The  latter  approximation  rests  on  the  supposition  that  the  contributions  in  (8)  from  terms 
involving  higher  derivatives  of  f(x)  are  small.  As  we  shall  see  later,  examination  of  the  result 
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obtained  in  this  way  shows  the  error  to  be  small  except  at  points  within  ar  of  the  ends  of  the 
cylinder.  This  approximation  is  justified  by  recognizing  that  we  have  neglected  the  terms 


g3r3  -i-  3 arx2 
3 


(13) 


This  series  for  f(x)  as  defined  in  (11)  is  clearly  convergent  when  x/1  is  small.  The  approximation 
adopted  is  acceptable  provided  the  ratio  of  successive  terms  is  sufficiently  small.  In  particular, 
it  is  necessary  that  the  ratio  of  the  first  neglected  term  is  small  compared  with  that  retained.  This 
ratio  is 


t*  a2r2  +  *2 
3  (/2  -*2)2 


(14) 


and  is  <  0.2  provided  that  x  <  0.75.  Similar  conclusions  follow  for  the  ratios  of  further  terms. 
Thus,  the  approximation  is  valid  everywhere  save  at  points  close  to  the  ends  of  the  whiskers. 
From  (5)  and  (7)  we  then  have 


ex  =  \ib2Par  ,df 
b  n(l -v)kT  dx 

whence 

F  -  [,(1~’1  M.  (f»  -  **)]■* 
\ib2  DAr 


and 


(15) 


(16) 


a(x)  * 


,  \iaekTr  m 
7t(l-v)6D 


x 

(l2  -  x2)* 


(17) 


cr(x)  is  a  shear  stress  which  acts  over  the  lateral  surface  of  the  cylinder  and  which  produces  a 
total  force 


2nrfla(x)dx 

Jo 


(18) 


We  equate  this  force  to  that  derived  from  a  uniform  normal  stress  which  acts  parallel  to  the  axis 
of  the  cylinder.  This  stress  is  given  by 
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Wxb  -  ILj*  <i9) 

7t(l  -v)b2Dr 

Then  on  the  basis  that  the  aspect  ratio  (1/r)  of  the  cylinders  is  large  we  suppose  that  the  normal 
force  in  the  cylinders  associated  with  this  stress  must  be  balanced  by  one  which  operates  in  the 
matrix  and  which,  to  a  first  approximation  is  equivalent  to  a  stress 

At  =  JL[  ]W  (20) 

1  ~c  it(l-v)b2Dr 


b.  Comparison  with  experiment 

To  examine  the  effect  of  this  stress  on  the  observed  strain  rate  we  recognize  that 

€  =  Kon  <21> 

where  a  is  the  operative  stress.  In  the  presence  of  fibers  the  operative  stress  becomes: 

x  -  Ax  =  x  -  ae 

On  substitution  of  best  estimates  of  the  values  appropriate  for  MoSi2;  y.  =  lOOGPa;  T  =  HOOK; 
1  =  30/im;  r  =3/im;  b=  0.3nm  and  with  20%  of  SiC  we  get,  almost  independent  of  e,  Ar  = 
5MPa.  This  value,  in  view  of  the  uncertainties  in  the  values  of  the  parameters,  is  in  excellent 
agreement  with  that  found  at  NRL  by  Sadananda13,  namely  10  MPa. 

c.  Implications 

★  The  analytic  theory  agrees  with  experiment. 

d.  Plans 

★  Numerical  computations  will  be  performed  to  check  the  influence  of  the  approximations  in  the 
analytic  theory. 

★  Comparison  with  experiment  will  be  extended  to  a  wider  range  of  environments. 
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Task  120.  Macroscopic  stochastic  model  for  creep 

The  methods  of  Boltzmann’s  classical  statistical  mechanics  have  been  employed  to  calculate  the 
yield  strength  in  the  absence  of  thermal  activation  in  the  case  where  hardening  is  due  to 
randomly  dispersed  spherical  particles. 

a.  Theory 

We  seek  to  determine  the  yield  strength  at  OK  of  a  material  hardened  by  the  presence  of 
spherical  particles  which  act  as  simple  blocks  to  the  motion  of  dislocations.  The  analysis  is  based 
on  a  consideration  of  the  stability  of  a  sequence  of  lengths  of  a  single  dislocation  which  has 
impinged  on  spherical  particles  in  the  process  of  traversing  a  slip  plane. 


We  assume  that  when  acted  on  by  a  uniform  stress  o,  dislocations  are  bent  into  arcs  of  circles 
of  radius  R,  where 


o  =  (23) 

R 

fi  is  the  shear  modulus  and  b  the  magnitude  of  the  Burgers  vector;  that  impenetrable  particles 
are  randomly  distributed  and  provide  restraining  forces  of  a  magnitude  which  is  related  through 
an  angle  a(r)  to  the  size  r  of  the  circle  of  intersection  of  particle  and  slip  plane  and  which  is 
specified  by 


F  =  pfc2sin  a (r)  @4) 

and  that  the  yield  stress  is  the  largest  for  which  it  is  possible  to  have  a  stable  dislocation  join  an 
infinite  number  of  particles  along  a  line  which  is  essentially  straight  and  hence  does  not  cross 
itself. 


Towards  the  determination  of  this  stress  we 
shall  adapt  the  theoretical  analysis  given  by 
Hanson  and  Morris14  and  by  Labusch15  for  the 
yield  strength  due  to  such  a  dispersion  of 
point  obstacles  all  of  which  provide  the  same 
maximum  restraining  force  of 

F  =  pt2sina  -  apt2  @5) 

This  analysis  rests  on  the  rolling  circle 
technique  employed  by  Foreman  and  Makin16 
in  their  numerical  investigation  of  this 
hardening.  Thus,  it  was  argued  (see  Figure  5) 


Figure  5  The  rolling  circle  construction 
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that  a  dislocation  stable  at  a  sequence  of 
positions  k-2,  k-1  will  also  be  stable  at  the 
next  position  k  provided  the  particle  associated 
with  the  k+ 1th  position  is  centered  anywhere 
in  the  region  traversed  by  a  rolling  circle  as 
shown  hatched;  this  region  is  bounded  by  an 
angle  0  ^  2a.  A  particular  characteristic  of 
such  sequences  of  dislocation  lengths  is  that 
they  must  not  cross  themselves. 


Here,  we  shall  adapt  the  treatment  of  Hanson 

and  Morris14  to  the  case  where  the  strength  of  a  particular  obstacle  depends  on  the  radius  of  the 
circle  of  intersection  between  the  spherical  particle  concerned  and  the  slip  plane.  To  proceed, 
we  note  first  from  Figure  5,  in  the  case  of  point  obstacles  that,  following  the  path  of  a  rolling 
circle,  the  angle  turned  between  successive  tangents,  tk,  is  equal  to  the  angle  <t>  which  gives  the 
direction  of  the  trailing  circle  edge  at  the  point  of  rotation  (see  Figure  6).  It  is  then  readily  seen 
that  in  the  case  of  an  obstacle  with  radius  r,  the  angle  turned  becomes 


4  =  <J>(1+-C)  =  r<j)  (26) 

R 

Then,  if  the  number  of  points  is  N  and  if  the  angle  4>  occurs  n(  times  we  have 

N  =  f>(.  (27) 

and  the  total  angle  turned  is 

®  =  E  nfii  (28) 

i-i 


If  N  is  indefinitely  large  but  <£  finite,  the  average  value  of  4>{  is  zero.  In  this  case  the  whole 
dislocation  line  is  either  straight  or  contains  equal  numbers  of  complete  clockwise  and 
anticlockwise  rotations.  The  latter  configuration  is  a  situation  of  low  probability  which  we  can 
neglect. 

The  procedure  for  the  determination  of  the  thermodynamic  probability  is  akin  to  that  employed 
to  determine  the  Maxwell-Boltzmann  distribution  of  kinetic  energies  in  an  ideal  monatomic  gas. 
Thus,  if  there  are  n  partitions  of  i  at  each  of  N  points  the  thermodynamic  probability  that  the 
ith  is  employed  n,  times  is: 
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P  = 


M  IlA ' 


1-1 


n 


n»i 


(29) 


where  p;  is  the  occupational  probability  of  the  ith  region.  To  determine  p,  we  divide  the  hatched 
area  A  into  n  regions  of  area  5 A*  (i  =  1,  n),  each  characterized  by  a  value  of  <f>,  say  and  see 
that  when  the  concentration  is  c  the  probability  that  a  point  k  lies  in  a  particular  area  5A;  is 

p{  =  1  -  exp(-C*A)  -  C6 A,  (30) 


provided  the  area  5A4  is  sufficiently  small.  On  the  basis  that  the  most  probable  configuration  of 
a  system  having  a  large  number  of  parts  is  that  which  actually  occurs,  to  determine  this 
configuration  we  have  only  to  maximize  P  subject  to  the  constraints  offered  by  (27)  and  (28). 
To  proceed  we  take  advantage  of  the  facts  that  N  and  p  are  both  large  to  employ  Stirling’s 
approximation  and  then 

n  n 

filnP  =  £  firt*  +  bv7/  _  £  hmi  (31) 

i-l  i-l 

while  (27)  and  (28)  become 


n 


E 


fin.  =  0 


(32) 


£  fifljfy  =  0 

Introducing  the  Lagrangian  multipliers  X  and  y  we  have 

lnp,  -  \nni  +  y  +  =  0 


(33) 


(34) 


irt*‘ 


i-l 


(35) 


then  from  (28)  and  (33)  we  have 
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N  ^ipi  =  <bS  (36) 

i-l 

£  *,  «*<•>  =  s  (37) 

i-l 

and  on  substituting  for  n;  in  the  logarithm  of  P  (31)  we  find,  after  some  manipulation 


P  =  sNe'k *  (38) 

Then  invoking  (30)  and  (35)  and  proceeding  to  the  limit  of  small  quantities  &4,  (38)  becomes 

P  =  [cfexidAf  (39) 


We  have  now  to  determine  s  and  X.  Here,  where  A  is  the  area  indicated  in  fig.  (1)  we  have 

2R  T\2a-tia~\lf2X)) 

S  -  cf  (/  +  r)dl  f  ek*d<X>  (4°) 

0  H&f'ilflK)) 


The  determination  of  X  is  assisted  by  the  realization  that 

dlnS  =  IdS  =  y  ^ 

dX  SdX  ft  S  "  N 


0  as  N  —  «• 


(41) 


and  X  must  be  such  that 


=  0  (42) 

dX 

This  question  can  be  resolved  only  by  numerical  means.  In  the  special  and  unrealistic  case  in 
which  a  <  1  we  find,  independent  of  the  value  of  T,  that  X  =  1.41/a. 

b.  Comparison  with  experiment 

*  In  comparison  with  the  case  treated  previously'6  there  are  additional  complications:  the 
obstacles  have  finite  radius  and;  because  the  radii  of  individual  obstacles  differ  one  from 
another,  they  have  variable  strengths.  Bacon,  Kocks  and  Scattergood17  have  shown  through 
numerical  studies  that  the  Orowan  stress,  r  for  a  linear  array  of  impenetrable  circular  blocks 
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to  dislocation  motion  having  radius  R  and  inter-center  spacing  L  is  given  by 

t  =  M  Jg*  (43) 

4n  L  -2R 

where  X  is  the  harmonic  mean  of  L  and  2R  (2).  Thus,  the  specific  strength  of  these  obstacles 
is  variable  and  denoted  by  the  term  lnX/2x.  Here,  the  obstacles  are  not  distributed  along  a 
straight  line  nor  regularly  spaced.  Additionally  their  strengths  are  variable.  This  additional  work 
must  be  done  before  a  comparison  with  experiment  can  be  performed, 

c.  Implications 

★  That  the  zero-temperature  is  analytically  tractable  means  that  we  have  a  good  measure  against 
which  to  check  more  general  numerical  models. 

d.  Plans 

*  The  results  will  be  used,  for  checking  purposes,  as  the  zero-temperature  limit  of  the  finite- 
temperature  model  planned  in  task  111. 
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Task  200.  Lifetimes  of  piezoelectric  ferroelectrics 

The  utility  of  PZT  as  a  ’smart  material’  stems  from  its  high  dielectric  constant,  which  gives 
large  changes  in  polarization  for  quite  small  electric  fields.  Unfortunately,  this  property  is  linked 
to  ferroelasticity,  which  can  be  destructive  of  the  mechanical  stability.  It  is  this  consequence, 
which  leads  to  premature  mechanical  failure  of  sensors  subjected  to  cyclic  loading,  which  forms 
the  critical  issue. 

a.  Theory 

In  order  to  determine  what  can  be  done  to  address  the  critical  issue,  it  is  instructive  to  review 
the  physics  of  ferroelectrics.  The  properties  of  PZT  are  primarily  a  consequence  of  a  phase 
transition  in  which  the  crystal  structure  changes  from  cubic  above  the  Curie  temperature  (about 
K)  to  body-centered  tetragonal,  with  the  c  axis  slightly  longer  than  the  a  axes  (0.5  to  3%, 
depending  on  composition).  The  cell-centering  Zr/Ti  atoms  are  positioned  slightly  off-center 
along  the  c-direction,  so  that  each  cell  has  a  small,  permanent  elecii.c  dipole  moment.  In  the  as- 
transfoimed  state  the  dipoles  are  arranged  in  small  domains,  with  a  volume  fraction  f  =  1/3 
oriented  along  each  of  the  three  crystal  axes,  so  that  there  is  no  net  polarization.  The  elastic 
energy  is 


U„  ■  -Vsc)2  (1*30(1  -fl  (44) 

4 

where  is  the  transformation  strain  and  f  is  the  volume  fraction  of  cells  oriented  with  c-axis 
parallel  to  a  particular  crystal  axis.  It  is  clear  that  the  elastic  energy  is  a  maximum  for  f  =  1/3. 
The  electric  energy  (the  interaction  energy  of  the  dipoles)  will  be  a  minimum  for  the  same  value 
of  f.  To  produce  a  useful  material,  the  material  must  be  polarized,  or  ’poled’;  that  is,  exposed 
for  a  short  time  (10-60  minutes)  to  a  large  electric  field  (of  order  1-10  MV/m)  at  a  temperature 
of  about  450K,  somewhat  below  the  Curie  temperature  but  above  the  service  temperature.  This 
causes  a  preferential  orientation  of  the  cell  c-axes  along  the  electric  field  axis,  in  a  sense  such 
that  the  dipoles  lie  as  closely  antiparallel  as  is  possible  to  the  applied  field.  This  decreases  the 
elastic  energy,  as  can  be  seen  from  (44)  if  f  is  taken  to  be  the  volume  fraction  of  dipoles  lined 
up  antiparallel  to  the  field  axis.  The  internal  electric  energy  is  increased,  but  the  total  electric 
energy  decreases  because  of  the  negative  energy  of  the  antiparallel  dipoles  in  the  applied  field. 
The  material  is  then  cooled  and  the  electric  field  removed,  leaving  a  low  elastic  energy,  a  high 
electric  energy  and  a  spontaneous  polarization  P,.  The  poled  state  is  never  maximally  polarized. 
The  best  that  can  be  done  is  to  rotate  the  dipoles  within  each  grain  to  the  crystal  axis  most 
closely  parallel  to  the  field  direction,  which  may  leave  a  dipole-field  angle  of  up  to  45°.  Even 
this  optimum  polarization  cannot  be  realized,  because  the  dipole  rotation  causes  grain  boundary 
stresses  to  develop  and  increase  to  the  point  that  further  rotation  is  suppressed.  Therefore  the 
poled  grains  consist  of  domains  which  are  differently  oriented,  but  with  a  large  net  dipole 
moment  in  the  direction  of  the  applied  field.  That  the  unpoled  and  poled  states  are  energetically 
metastable  indicates  that  there  is  an  energy  barrier  associated  with  the  change;  otherwise  a 
transition  to  the  ground  state  would  occur  spontaneously.  Therefore  changes  in  the  polarization 
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take  place  by  one  or  more  thermally  activated  processes,  for  example  domain  nucleation  and 
growth  or  domain  wall  motion.  The  magnitude  of  P,  after  poling  is  of  order  0.1-  0.5  coulombs 
m'2.  The  electrostatic  energy  density  of  a  system  of  parallel  dipoles  in  SI  units  is 


Po  y  3cos20i  -  1 
8*Va  «•  r;3 


(45) 


where  p0  is  the  dipole  size,  k0  is  the  permittivity  of  free  space,  v,  is  the  volume  of  the  unit  cell. 
The  sum  is  performed  for  a  single  dipole  and  is  taken  over  all  neighboring  dipoles,  with  r;  the 
magnitude  of  the  vector  joining  the  origin  to  the  ith  neighbor  and  6{  the  angle  between  this  vector 
and  the  dipole  axis.  For  a  cubic  structure  the  sum  in  (45)  vanishes  by  symmetry;  therefore  the 
magnitude  of  the  electrostatic  energy  depends  on  the  transformation  strain  5e  (in  fact  the 
dependence  varies  as  5e2,  the  same  as  the  elastic  energy).  For  8e  of  0.01,  the  sum  in  (45)  has 
a  value  of  about  0.  lv,. 

For  a  transformation  strain  of  0.01  and  a  shear  modulus  of  lOOGPa,  the  elastic  energy  of  the 
unpoled  state  is  about  3MJ/m3.  In  a  completely  poled  state,  with  P,  =  0.1  -  0.5C/m2,  the 
electrostatic  energy  varies  from  5  to  125MJ/m3.  Although  these  are  only  approximate  quantities, 
the  indication  is  that  the  electrostatic  energy  stored  during  the  poling  process  is  substantially 
larger  than  the  lost  elastic  energy.  This  means  that  the  poled  material  exists  in  a  relatively  high- 
energy,  metastable  state. 

Now  attention  will  be  turned  to  the  question  of  premature  failure  under  cyclic  deformation, 
which  occurs  under  conditions  of  mechanical  loading18  for  quite  moderate  stress  levels.  The 
propensity  toward  mechanical  instability  is  not  surprising  in  view  of  the  magnitude  of  the  stored 
electrostatic  energy,  but  the  mechanism  has  yet  to  be  determined.  The  specific  experiments 
cited18  involve  cyclic  bending  such  that  the  tensile  and  compressive  stresses  developed  act  normal 
to  the  poling  axis.  The  larger  part  of  the  specimen,  therefore,  consists  of  domains  in  which  the 
polarization,  together  with  the  c-axis,  is  oriented  more  normal  than  parallel  to  the  stress.  In  the 
tensile  part  of  the  cycle,  domain  walls  will  migrate  to  favor  domains  oriented  more  parallel  to 
the  stress  direction  and  new  domains  will  nucleate  and  grow.  This  effect  is  analogous  to  the 
poling  process  in  that  the  field-coupled  elastic  energy  is  increased,  but  is  compensated  by  a 
negative  work  done  by  the  applied  stress.  At  the  same  time  the  electrostatic  energy  is  reduced. 
When  the  tensile  stress  is  released,  the  stored  elastic  energy  becomes  manifest  as  internal 
stresses,  principally  at  the  grain  boundaries.  Because  the  electrostatic  energy  in  the  poled  state 
is  so  large,  the  elastic  energy  developed  can  be  much  larger  than  in  the  original,  unpoled  state 
and  failure,  most  likely  by  intergranular  fracture,  is  inevitable. 


b.  Comparison  with  experiment 

The  arguments  above  explain  the  basic  physical  reasons  for  failure,  but  not  the  specific 
mechanics.  In  the  compressive  cycle,  the  effects  discussed  above  will  be  reversed,  which  negates 
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the  simple  rationalization.  Therefore  two  observations  remain  to  be  explained. 

1 .  According  to  the  mechanism  above,  a  simple  monotonic  tensile  stress,  applied  for  a  sufficient 
time  and  then  released,  would  be  adequate  for  failure  to  occur  by  thermally  activated  creep  of 
the  domain  structure.  There  are  no  explicit  experimental  measurements  to  confirm  or  refute  this 
hypothesis.  However,  it  is  reasonable  to  expect  that  the  effect  would  have  been  noticed  and 
reported,  if  it  occurred. 

Domain  walls  are  susceptible  to  pinning  by  defects19;  these  may  be  vacancies,  substitutional 
inhomogeneities  or,  of  course,  grain  boundaries.  It  is  reasonable  to  expect  that  the  response  to 
a  monotonic  stress  will  be  decay  into  a  state  in  which  all  walls  are  pinned,  so  that  further 
development  of  the  substructure  is  inhibited.  In  fact  the  experimental  results  from  NIST18  show 
that  the  fracture  stress  at  room  temperature  in  monotonic  bending  is  about  140MPa,  or  about  1.4 
x  10'V  This  is  a  relatively  low  stress  and  may  reflect  a  catastrophic  release  of  energy  which 
is  likely  to  occur  once  a  crcck  is  nucleated. 

2.  Cyclic  stressing,  according  to  the  same  mechanism,  should  be  completely  reversible. 
However,  experiment18  shows  that  under  cyclic  stressing  at  70%  of  the  fracture  stress,  the 
lifetime  is  in  excess  of  107  cycles.  At  stress  amplitudes  from  75  %  to  95  %  of  the  fracture  stress, 
the  lifetime  decreases  sharply  from  105  to  102  cycles  (although  the  spread  is  large).  Clearly,  not 
only  is  there  is  a  threshold  stress,  but  there  is  also  a  cyclic  irreversibility  for  stress  levels  in 
excess  of  the  threshold. 

The  presence  of  a  threshold  stress  is  compatible  with  thermally  activated  domain  wall  motion. 
There  are  several  reasons  for  irreversibility.  The  most  obvious  is  that  the  processes  of  domain 
annihilation  and  nucleation  are  not  reversible.  Another  stems  from  the  pinning  of  domain  walls, 
as  discussed  above.  Once  a  wall  segment  is  pinned,  it  cannot  easily  retrace  its  propagation  path. 
Therefore  wall  motion  in  the  depoling  sense  will  not  easily  reverse,  instead  proceeding  in  the 
depoling  sense  by  a  ratchet  mechanism.  The  depoling  is  associated  with  an  increased  elastic 
energy  density,  which  will  be  manifest  as  steadily  increasing  internal  stresses  each  time  the 
applied  cyclic  stress  is  relaxed.  A  third  reason  for  irreversibility  results  from  the  binding  of 
dislocations  to  the  domain  walls.  With  a  misfit  of  1%  at  the  interface  between  parallel  and 
normal  domains,  a  wall  of  length  just  a  hundred  lattice  constants  will  accumulate  a  total  misfit 
of  a  complete  lattice  constant;  the  energy  of  the  wall  can  then  be  reduced  substantially  by 
forming  a  misfit  dislocation,  that  is  an  edge  dislocation  with  Burgers  vector  in  the  plane  of  the 
wall  and  with  the  extra  half-plane  lying  always  in  the  normal  domain.  Then  wall  motion  is 
limited  by  the  rate  at  which  the  dislocations  can  climb.  In  the  tensile  part  of  the  cycle,  the 
dislocations  must  climb  by  emitting  interstitials  (or  absorbing  vacancies),  while  in  the 
compression  phase,  vacancies  are  emitted.  Any  asymmetry  in  these  two  processes  will  be 
manifest  as  an  irreversibility.  The  advantage  of  this  climb  mechanism  is  that  the  predicted 
irreversibility  is  cumulative,  a  property  not  shared  by  the  other  irreversible  processes  discussed 
above.  Finally,  it  is  also  important  to  look  at  the  piezoelectric  properties  of  PZT,  even  though 
the  case  in  question  involves  mechanical  loading  at  constant  (=0)  electric  field.  Consider  the 
simple  case  in  which  the  tensile/compressive  stress  an  lies  along  the  x,  axis  and  the  poling  axis 
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is  x3.  Then  there  is  a  piezoelectric  coupling 

=  ^311°  n  (46) 

where  Aps  is  the  polarization  increment  and  dijk  is  the  piezoelectric  tensor.  The  coefficient  d3n 
is  negative,  so  that  the  dipoles  lined  up  antiparallel  to  the  poling  axis  are  increased  in  magnitude. 
The  dipoles  normal  to  the  poling  axis,  however,  are  unchanged.  Therefore  the  electrostatic 
energy  density  becomes  larger  in  the  parallel  (to  the  poling  axis)  domains  than  in  the  normal 
domains.  This  means  that  domain  wall  motion  of  such  a  sense  as  to  increase  the  size  of  the 
normal  domains  (as  is  the  case  in  tension)  is  accompanied  by  a  small  decrease  in  electrostatic 
energy  and  will  be  enhanced.  In  compression  the  energy  in  the  parallel  dipoles  becomes  smaller. 
Since  the  stress  is  reversed,  the  wall  motion  is  reversed  and  is  similarly  enhanced.  Thus  the 
piezoelectric  effect  is  completely  reversible. 

c.  Implications 

★  We  conclude  that  the  irreversibility  which  leads  to  failure  of  PZT  in  cyclic  stressing  is 
connected  with  ferroelastic  domain  wall  motion  and  is  not  influenced  by  the  piezoelectric 
character.  If  the  problem  is  to  be  corrected,  it  is  essential  to  determine  the  mechanism  which 
controls  domain  wall  mobility.  To  this  purpose,  both  theory  and  experiment  are  required. 

d.  Plans 

★  The  mechanisms  of  ferroelastic  domain  wall  motion  will  be  modelled,  with  the  aim  of 
determining  critical  parameters  which  can  be  derived  from  or  compared  with  experiment. 
Specific  questions  to  be  answered  are: 

1.  Is  the  wall  motion  governed  by  the  two-dimensional  analog  of  kink  pair  nucleation  and 
propagation,  or  by  diffusive  dragging  of  wall  defects  and  dislocations? 

2.  How  are  these  mechanisms  influenced  by  an  electric  field  of  sufficient  magnitude  to  stabilize 
the  poled  state? 

3.  How  do  the  lifetimes  of  poled  and  unpoled  specimens  compare? 

★  The  theoretical  work  will  be  closely  coordinated  with  experimental  work  at  NIST  and  NRL. 
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Figure  Captions 

Figure  1  Definition  of  the  geometric  parameters  D,  L  and  0. 

Figure  2  An  Orowan  plot  of  strength  for  dispersion  hardened  materials. 
Figure  3  Plot  of  strength  following  the  unified  approach  of  equation  (2). 
Figure  4  The  distribution  of  breaking  angles. 

Figure  5  The  rolling  circle  construction. 

Figure  6  Relationship  between  and  t^. 
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